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2^ ■ Abstract 



We present a general purpose Monte Carlo program for the calculation of any 
infrared safe observable in Wj and Z7 production at hadron colliders at next- 
to-leading order in a s . We treat the leptonic decays of the W and Z-boson in 



the narrow-width approximation, but retain all spin information via decay angle 
correlations. The effect of anomalous triple gauge boson couplings is investigated 
and we give the analytical expressions for the corresponding amplitudes. Fur- 
thermore, we propose a way to study the effect of anomalous couplings without 
introducing the ambiguity of form factors. 



1 Introduction 



The production of W'y and Z*y in hadronic collisions has been studied extensively since 
the Born cross sections have been computed |l|, 0. In particular, these processes allow 
to study the triple gauge boson couplings WW'j, ZZj and Z77. The study of these 
couplings is mainly motivated by the hope that some new physics may modify them. If 
the new physics occurs at an energy scale well above that being probed experimentally, 
it is possible to integrate it out. The result is an effective theory which might result in 
non standard triple gauge boson couplings. 

Both collaborations at the Tevatron have studied the production of Wj || and 
Z7 0] pairs. The bounds on the anomalous couplings obtained at the Tevatron tend 
to be less constraining than those obtained at LEP |J. However it has to be kept in 
mind that these analyses are complementary. At the hadron colliders a whole range 
in the center of mass energy is tested, whereas at LEP the center of mass energy is 
fixed by the collider. Furthermore, with Run II the expected number of events at the 
Tevatron increases substantially. Assuming a data sample of 2 fb _1 , more than 3000 
W'y — > £wy events and 700 Z7 — > ^£7 events are expected for each experiment ||. Of 
course, the expected number of events is even bigger for the LHC. 

Anomalous triple gauge boson couplings lead to deviations from Standard Model 
predictions. Obviously, observables for which these deviations are enhanced offer better 
chances to find new physics or get tighter constraints on anomalous couplings. There 
are basically two classes of such observables. Either we consider observables which 
are strongly suppressed in the Standard Model or observables with large transverse 
momentum (or center of mass energy). In both cases, the inclusion of next-to- leading 
order (NLO) QCD corrections is mandatory. 

A prominent example of an observable that is suppressed in the Standard Model 
is the so called radiation zero for W^y production. At leading order (LO) there exist 
some kinematic configurations for which the amplitude vanishes [TJ. This is manifest 
in some observables as a dip in the rapidity distributions. Since anomalous coupling 
contributions fill in the dips, there seemed to be excellent prospects to obtain accurate 
limits for them from experimental data. Unfortunately, next-to-leading order QCD 
corrections strongly affect the LO analysis. They have the same effect as the anomalous 
coupling contributions. The dips are filled in, making the extraction of anomalous 
couplings quite more difficult. 

For processes with large transverse momentum or center of mass energy, the NLO 
corrections are particularly large. This is due to the fact that the cross sections in these 
cases get large contributions from gluon induced partonic subprocesses, which only 
enter in a next-to-leading order description of the cross section. Thus, even though the 
anomalous contributions are enhanced in these regions, a calculation at NLO in a s is 
required to reliably exclude (or establish) physics beyond the Standard Model. 

The relevance of NLO corrections was first shown for the production of real (spin- 
summed) W and Z bosons with Standard Model couplings and without considering 
lepton decays and spin correlations |7], |]. These calculations were later extended, in 
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order to include the leptonic decays and anomalous couplings [§, [K| |Tl|]. However, the 
full one-loop amplitudes including leptonic decays became available only very recently 
12fl . Therefore, refs. |9|, |ll|] included decay correlations everywhere except for the 



finite part of the virtual contributions. 

In this paper we present order a s results for the production of Wj and Zj in 
hadronic collisions, including the full leptonic correlations. We work in the narrow- 
width approximation, where only 'single-resonant' Feynman diagrams have to be con- 
sidered. The simplicity of the helicity method allows to take into account anomalous 
couplings as well and present for the first time analytical expressions for the corre- 
sponding amplitudes. 

For the case of WW, ZZ and WZ production at hadron colliders, some results 
beyond the narrow-width approximation are known. The narrow-width approximation 
requires only the calculation of 'doubly-resonant' Feynman diagrams. However, for 
these processes, also the amplitudes including 'single- resonant' diagrams have been 



computed and implemented into a Monte Carlo program [[13 



To perform the phase space integration we use the subtraction method discussed 
in ref. fl4| |. This allows for a straightforward implementation of the one- loop qq' — * 



V^i — > ££'7 and bremsstrahlung qq' — > gVj — > gii'j amplitudes, presented in ref. (12 
(V E {Z,W}). The constructed Monte Carlo code allows the computation of any 
infrared-safe observable^. 

A brief overview of the calculation is given in section 2, were we summarize the input 
parameters used and the cuts implemented to obtain our phenomenological results. In 
section 3 we first present some benchmark cross section numbers for both W'j and Z'j 
production at the LHC and study the typical scale dependence of some observables 
at NLO in the Standard Model. Since many distributions have been studied in the 
past, we refrain from doing a detailed analysis. However, as soon as more precise data 
becomes available such an analysis can easily be done. 

In section 4 we concentrate on anomalous triple gauge boson couplings. We describe 
the parameterization of the triple gauge boson vertex in terms of anomalous coupling 
parameters and search for the kinematical region where its effect is amplified, namely 
at large transverse momentum for both photon and leptons. We also analyze the 
possibility of seeing the effect of approximate radiation zeros in the W^ process, i.e. by 
looking for 'dips' in rapidity distributions. In order to avoid the arbitrariness introduced 
by form factors, we propose to analyze the anomalous couplings as a function of the 
squared partonic center of mass energy s. This has been suggested previously for 
the Z7 case ||15|| , where such an analysis is straightforward. We extend this idea to 
the Wj production. This case is more involved, since a complete reconstruction of 
s is impossible, due to the appearance of a non observed neutrino in the W decay. 
Particularly, we present an observable quantity which is highly correlated to s and, 
therefore, allows such an analysis even for W^ production. Finally, in section 5 we 
give our conclusions and in the appendix we present analytical expressions for the 



1 The corresponding Fortran codes are available upon request 
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amplitudes relevant for the processes under consideration. 



2 Formalism 

The helicity amplitudes needed for the computation of the NLO corrections to W^f and 
Z7 production in the Standard Model were presented in ref. [T2|. The amplitudes rele- 



vant for the inclusion of anomalous couplings are presented in the appendix. In order to 
cancel analytically the soft and collinear singularities coming from the bremsstrahlung 
and one loop parts, we have used the version of the subtraction method presented in 
ref. [ nj . The amplitudes are therefore implemented into a numerical Monte Carlo style 



program, which allows to calculate any infrared-safe physical quantity with arbitrary 
cuts. 

Obviously, the Monte Carlo program can be used for the Tevatron and the LHC. 
However, in this paper we will mainly concentrate on results for the LHC collider, which 
corresponds to pp scattering at ^fs = 14 TeV. Unless otherwise stated, the results are 
obtained using the following cuts: we make a transverse momentum cut of p l T > 25 GeV 
for the charged leptons and the rapidity is limited to \r]\ < 2.4 for all detected particles. 
The photon transverse momentum cut is pj. > 50(100) GeV for W^f (Z7) production. 
For the Wj case we require a minimum missing transverse momentum carried by 
the neutrinos p™ lss > 50 GeV. Additionally, charged leptons and the photons must 
be separated in the rapidity-azimuthal angle by AR^ = yj (?7 7 — r]i) 2 + (0 7 — (pi) 2 > 
0.7. Moreover, since our calculation is done in the narrow- width approximation and, 
therefore, ignores the radiation of photons from the final state leptons, we apply an 
additional cut to suppress the contribution from the off-resonant diagrams. For that 
purpose, we require the transverse mass M? > 90 GeV for W'y production and the 
invariant mass of the ££7 system M« 7 > 100 GeV for the Z'j case. 

Finally, photons can also be significantly produced at LHC from the fragmentation 
of a final state partonP| Unfortunately, fragmentation functions of partons into photons 
are not very well determined and the NLO calculation for such a contribution is not 
available yet. In principle, a full NLO calculation should include it however, since only 
the sum of the 'direct' plus 'fragmentation' components is physically well defined at 
NLO (only in the sum all collinear singularities cancel out). In order to circumvent 
this problem, we include the LO component of the fragmentation part but using NLO 
fragmentation distributions, where we can factorize the final state 57 collinear singu- 
larities. Since the 'fragmentation' component can be further suppressed implementing 
certain cuts (see below) the lack of its NLO calculation is not expected to affect the 
final result beyond the few percent level. 

The fragmentation contribution constitutes a background to the search of anoma- 
lous couplings, since it does not involve any triple gauge boson coupling. Fortunately, 
there is a way to suppress its contribution by requiring the photons to be isolated from 
hadrons. In this paper we require the transverse hadronic momentum in a cone of 



2 This contribution is also known in the literature as 'bremsstrahlung' 
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size Ro = 0.7 around the photon to be smaller than a small fraction of the transverse 
momentum of the photon 

E PT d <0-15^ (1) 

AR<R 

This completes the definition of our 'standard' cuts. 

When indicated, we also apply a jet-veto, which means that we reject any event 
where a jet of > 50 GeV and |?j et | < 2.5 is observed. 

In our results we do not include the branching ratios of the vector boson into 
leptons. They can be taken into account by simply multiplying our final results with 
the corresponding branching ratio. For both the LO and NLO results we use the latest 
(corrected) set of parton distributions of MRST(corOl) |L6| and the two loop expression 
(with rif — 5 for the typical scales of these processes) for the strong coupling constant 
with Aj^g(ny = 4) = 300 MeV which corresponds to a s (M z ) = 0.1175. For the 



fragmentation component we use the fragmentation functions from ref. [17]. 

Since we are particularly interested in the large px tail, which is more sensitive to 
the anomalous coupling contributions, we use (unless otherwise stated) 

/i 2 = £ = My + ^ [{Plf + (pi) 2 ] . (2) 

as the 'standard' scale for both the factorization and renormalization scales. 

Contributions from b and t quark initial states have been neglected and, con- 
sistently, the following values have been used for the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements in the case of production: \V u d\ = \V CS \ = 0.975 and 
\V US \ = \V cd \ = 0.222. 

The masses of the vector bosons have been set to Mz = 91.187 GeV and Mw = 
80.41 GeV. We do not include any QED or electroweak corrections but choose the 
coupling constants a and sin 2 #jy in the spirit of the 'improved Born approximation' 
]TR with sin 2 9w = 0.230. Notice that the observable is order a 2 ; within the same 



spirit we use the running a = a(Mz) = 1/128 for the coupling between the vector 
boson and the quarks (to effectively take into account the EW corrections) whereas we 
keep a = 1/137 for the photon coupling. It is worth noticing that this modification 
results in a 7% change in the normalization of the cross section with respect to the 
standard approach of using both running coupling constants. 



3 Standard Results 

We begin the presentation of our results with some total cross section numbers in 
Table [I], which can be useful for future checks and for an estimate of the number of 
events to be observed at the LHC The first three results were obtained by impos- 
ing only the cut on the transverse momentum of the photon, i.e. pj. > 50(100) GeV 
for W"y (Zj) production. Apparently, the NLO corrections as well as the fragmenta- 
tion contribution are very large. As discussed above, the relative importance of the 
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fragmentation contribution can be reduced substantially by applying the isolation cut 
prescription. This can be seen from the results for the total cross section obtained after 
the implementation of our standard cuts, which are also presented in Table [l|. Unfor- 
tunately, most of the previous publications on the subject [§, |], [10], |TTJ do not present 
cross sections numbers. Nevertheless, we have compared many of the plots shown in 
refs. m H |TTJ, specially for the case of real (spin-summed) W/Zj production f| 
(which is not affected by lepton correlations). Within the precision that can be reached 
in such a comparison, we found good agreement. 



a 


LO* 


Frag. 


NLO 


W+7 {pi > 50 GeV) 


4.79 


3.02 


13.89 


W"7 if T > 50 GeV) 


3.08 


3.55 


10.15 


Z-i (pi > 100 GeV) 


1.29 


0.412 


2.37 


W + ~i (std. cuts) 


0.436 


0.094 


1.71 


(std. cuts) 


0.310 


0.095 


1.20 


Z'-f (std. cuts) 


0.524 


0.041 


0.877 



Table 1: Cross sections in pb for pp collisions at = 14 TeV for fi = /i st . The statistical 
errors are ±1 within the last digit. LO* corresponds to the direct component only. 

In what follows we will estimate the theoretical uncertainty of our results by ana- 
lyzing the variation of various distributions when changing the scale /i by a factor of 
two in both directions l/2/i st < /i < 2/i st . Since many observables already have been 
studied in the past || ^ |K], [HI an d i n order to avoid the proliferation of plots, we 
refrain from presenting a detailed analysis here. We simply concentrate on a couple 
of typical examples in order to give a general picture and illustrate the importance of 
NLO corrections in both Wj and Zj production. 

In Figure [I] we show the scale dependence of the pr distribution of the photon 
in W + j production with the standard cuts applied (upper curves) and also with the 
additional requirement of a jet- veto (lower curves). As can be observed, the scale 
dependence is still large ( ± 10%) as long as only the standard cuts are applied. 
However, it is considerably reduced when the jet-veto is applied. The situation is 
similar to what has been observed in the case of WW production [p(J and is caused by 
the suppression of the contribution from the qg initial state appearing for the first time 
at NLO. Since this initial state dominates the cross section, the NLO result behaves, 
regarding the scale dependence, effectively like a LO one. 

In the inset plot we present the ratio between the NLO and LO results (with the 
standard scale), which remains larger than 3 and increases with the photon transverse 
momentum. This clearly shows that the LO calculation is not even sufficient for an 
understanding of the shape of the distribution, since the NLO effect goes beyond a 
simple change in the normalization. As is well known ||, the relevance of the NLO 
corrections for this process is mainly due to the breaking of the radiation zero appearing 
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at LO and to the large qg initial state parton luminosity at the LHC We also show 
the ratio of the NLO jet-veto and the LO result. As expected, this ratio is closer to 1, 
again due to the fact that most of the contributions coming from the new subprocesses 
appearing at NLO are suppressed by the jet-veto. 

It is worth mentioning that the scale dependence of the LO result turns out to be 
very small. This is an artificial effect and illustrates that a small scale dependence is 
by no means a guaranty for small NLO corrections. In fact, there is no renormalization 
scale dependence at all at LO. The only scale dependence comes from the factorization 
scale dependence of the parton distribution functions. Furthermore, we would like to 
mention that the situation concerning the scale dependence is slightly more favorable 
at the Tevatron. This is simply due to the fact that the gluon initiated process is less 
important. 




200 400 600 800 



P? (GeV) 

Figure 1: Scale dependence of a^ NL °^ without (upper curves) and with (lower 
curves) jet-veto. The scale has been varied according to ^(dashes) < (j, < 
2// s t(dots). The inset plot shows the ratio a^- N LO " > / a^ LO \ again without (solid) 
and with (dots) jet- veto. 

In Figure |2| we study the lepton correlation in the azimuthal angle for Z7 produc- 
tion, A(pu = \<p£- — Notice that this observable can be studied at NLO since we 
take fully into account the spin correlations between the leptons in our implementation 
of the one-loop corrections. The NLO corrections are rather sizeable and increase the 
cross section by 50% for small Acj) U . The region A(f> u > 2 (with our standard cuts) 
is kinematically forbidden unless a jet with a high transverse momentum is produced. 
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Therefore, the cross section vanishes at LO and it is strongly suppressed for the NLO 
calculation with jet-veto. In this region, the full NLO calculation is effectively only a 
LO calculation and its scale dependence becomes larger, as expected. 

Because there is no radiation zero appearing at LO for Z7 production, the NLO 
corrections are under better control in the kinematical region where the LO cross 
section does not vanish. Nevertheless, for large transverse momentum, the process with 
a qg initial state again dominates the NLO contribution and the corrections increase 
considerably [|K|. 




Figure 2: Scale dependence of a^ NL °^ without jet- veto (upper solid curves), 
a (NLO) w ith jet- veto (lower solid curves) and a^ L °^ (dotted curves). The scale 
has been varied according to ^ < /x < 2/i s t- 



Finally we mention that we also considered a more stringent photon isolation pre- 



scription, introduced by Frixione |21J. This prescription completely eliminates the 
fragmentation contribution. We have checked that the main features of all studied 
distributions remain unchanged when it is imposed. 



4 Sensitivity to Anomalous Couplings 

The study of triple gauge boson couplings is motivated by the hope that some physics 
beyond the Standard Model leads to a modification of these couplings which eventu- 
ally could be detected. In order to quantify the effects of the new physics an effective 
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Lagrangian is introduced which in principle contains all Lorentz invariant terms. The 
prefactor of these operators are the anomalous couplings. A general approach is im- 
practical, since it would lead to a proliferation of new couplings. Therefore, some 
additional constraints have to be imposed. 

The usual choice for W'-f production is to insist on electromagnetic gauge invariance 
and on C and P invariance. Also, only operators of dimension six or less are considered. 
This leads to a momentum-space vertex W~(q)Wt(q)'y ft (p) (where all momenta are 
outgoing p + q + q = ) which can be written as [G^, 120] 



7 t 



M 2 



r^ 7 (g, q,p) = q a 9^(2 + Ak 1 + A 7 ^-) - (2 + A/c 7 + A 



^ (1 + ^) + ^ V 



(3) 



where the overall coupling has been chosen to be — |e|. Note that in the Feynman rule 
for this vertex there is also a factor i that is conventionally not included in T Q/3/ \ In 
the Standard Model we have Ak 7 = A 7 = 0. 

For Z~f production, we consider operators up to dimension 8 (all of them C odd) 
resulting in ZZ^ and Z77 couplings. The non-standard Z a {q\)^p{q2)Z ll {p) momentum- 
space vertex is given by 



Ml V M z 



hle^q 2u - h^ e ^p a p v q 2 A (4) 



Ml 



where the overall coupling has been chosen to be \e\ (and e 0123 = +1). This ver- 
tex is absent altogether in the Standard Model. The non-standard -^(^^((^^(p) 
momentum-space vertex can be obtained from eq. (^) by setting q\ — ■> and replac- 
ing hf — > h] . Notice that the vertex differs from the one implemented in ref. [II 



by an overall factor i, which ensures the hermiticity of the corresponding effective La- 
grangianP] fl5| . Furthermore, the % factor modifies the interference pattern of anomalous 
coupling and Standard Model amplitudes: CP-violating h\,i contributions do not inter- 
fere with the Standard Model ones, whereas the CP-conserving ^4 do [[15]. Therefore, 
with the corrected vertex, there are contributions linear in h^^ and the cross section is 
generally not invariant anymore under a change of sign of h 3 ^. This must be considered 
in a future precise analysis of anomalous couplings from experimental data since the 
limits for /i 3 4 will not be symmetric at variance with present analyses. 

The anomalous couplings spoil the gauge cancellation in the high energy limit and, 
therefore, will lead to violation of unitarity for increasing partonic center of mass 

3 There is also a sign difference in the ^4 contributions coming from the different definition of 
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energy Vs. Usually, in an analysis of anomalous couplings from experimental data 
in hadronic collisions this problem is circumvented by supplementing the anomalous 
couplings, axe-, with form factors. A common choice for the form factor is 

where n has to be large enough to ensure unitarity and A is interpreted as the scale for 
new physics. Obviously, this procedure is rather ad hoc and introduces some arbitrari- 
ness ||23|| . Also, it is not really consistent with the effective theory approach. Increasing 
anomalous contributions would require the inclusion of even higher dimensional oper- 
ators. At the end of this section section we will address the question on how to avoid 
this arbitrariness in an analysis of anomalous couplings at hadron colliders. 



InB' 



Figure 3: Ratio of <tac an d ctsm at NLO for fj, = yu st . The anomalous couplings 
in o"ac have been chosen as Ak 7 = 0.08, A 7 = 0.02 and a form factor defined in 
eq. (|5]) with n = 2 and A = 2 TeV has been used. 



Anomalous couplings mainly affect the events with large s or large pr- Since the 
total cross section is dominated by low px events this is not a good observable to get 
tight constraints on anomalous couplings. A more promising possibility is to consider a 
double binned cross section. We therefore consider the total cross section binned in p T 
and p T for the process pp — > — > t + v~y at the LHC. In Figure we show the ratio 
of cr AC over <7sm where <tsm is the Standard Model cross section and <tac is the cross 
section obtained with An = 0.08 and A = 0.02 (both within the present experimental 
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limits from LEP || and Tevatron [|4]]) and a form factor as defined in eq. @ with 
n = 2 and A = 2 TeV. As expected, the ratio is large for the high bins, whereas it 
is very close to one for the low pr bins. We checked that the uncertainty coming from 
the scale variation is much smaller than the effect of the anomalous couplings for the 
high px bins. 

Another possibility to get a large effect due to anomalous couplings is to consider 
the approximate radiation zeros present in the Wj process [|l|]. At tree- level in the 
Standard Model, the Arj Wl = rjw — V-y distribution has a radiation zero. This dip is 
filled by next-to-leading order corrections and anomalous effects. In order to get an 



observable quantity, we do not consider Ar/wy but rather Ar] y g [^5j . This will wash out 
the dip as the rapidity of the lepton is not equal to the rapidity of the W. However, 
requiring the energy (or the transverse momentum) of the lepton to be large enough 
forces the lepton to follow closely the W direction. Also, applying a jet-veto reduces 
the effects of the next-to-leading order corrections. Thus, for larger p T even the next- 
to-leading order Arj^i distribution shows a clear dip in the Standard Model. This is 
illustrated in Figure |], where we show the Ar) 7 e distribution for the standard cuts and 
the additional cuts p T > 100 GeV and p f T > 200 GeV respectively. In all figures we 
apply our jet-veto. Also shown are the three distributions with anomalous couplings, 
which are chosen as for Figure || Clearly, the effect is dramatic for high energy leptons 
but of course, the disadvantage of applying such a cut is a big loss in statistics. 
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Figure 4: Ary 7 | distribution at NLO without (dashed curves) and with anoma- 
lous couplings Ak 7 = 0.08, A 7 = 0.02 (solid curves). In the first plot we applied 
only standard cuts and the jet-veto, in the second plot there is an additional 
cut p T > 100 GeV, and in the third plot p T > 200 GeV. 



We now turn to the question whether it is possible to avoid using form factors in the 
analysis of anomalous couplings at hadron colliders. This would bring these analyses 
more into line with those at e + e~ colliders. In order to do so one should analyze the 
data at fixed values of s, as it is done at LEP. This results in limits for the anomalous 
parameters which are a function of s. 
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Figure 5: Cross section at NLO for fj, = /i st with standard cuts and > 
200 GeV and p^ > 200 GeV. (a) no anomalous couplings; (b) hj = /if = 
0.01, hj = hf = 10 -4 and the usual dipole form factor with A = 2 TeV; 
(c) hj = h§ = 0.01, hj = h% = 10" 4 and no form factor; (d) hj = h$ = 
0.001, hj = hf = 10~ 5 and no form factor. The inset plot shows again (b) and 
(b) with the opposite sign in the anomalous couplings. 



Obviously, it is possible to do such analysis for the production of when both 
leptons are detected. Since the center of mass partonic energy can be reconstructed 
from the kinematics of the final state particles^ the cross section can be measured for 
different bins of fixed s ||15|| . As an example, we show in Figure |5| the cross section as a 
function of s. In order to enhance the effect of the anomalous couplings we do not only 
apply our standard cuts but we also require > 200 GeV and > 200 GeV. We show 
four curves: the curve (a) is the Standard Model result; curve (b) includes anomalous 
couplings in the standard way, that is with a form factor as defined in eq. (|Sp with n = 3 
(n = 4) for h\{hl) and A = 2 TeV and we have set hj = /if = 0.01, hj = hf = 10~ 4 ; 
curve (c) uses the same values for hf^" 1 but does not include any form factor; finally 
curve (d) does also not include any form factors but the anomalous couplings are smaller 

= /if = 0.001, hj = hf = 10~ 5 . For all of them we set h\ = h\ = 0. Of course, 
for large s the effects are much more dramatic if the form factor is omitted and at 
some point the corresponding curves would violate unitarity. This simply reflects the 
breakdown of the effective theory approach in this region. The idea behind Figure |5] is 
that such an analysis can be done for suitably defined bins in s. As a result, for each 



To simplify the discussion we assume that all final state particles are detected, including the jets. 



11 



bin, i.e. each value of s a bound on the anomalous couplings can be obtained. 

In the inset plot we present again curve (b) and the one corresponding to the 
parameters of (b) but with the opposite sign for values of the anomalous couplings 
hi = h§ = -0.01, hj = hf = — 10~ 4 . From there the effect of the interference 
between AC and SM results can be observed, i.e. the contribution of linear terms 
in h\ and h\ appearing in the squared amplitudes due to the correct treatment of 
the i factor in the ZjV vertex ||15|| . For this particular configuration, the interference 
effects are mostly relevant at V§ < 1 TeV and modify the cross section by more than 
10%. Clearly, they must be taken into account in a precise extraction of anomalous 
couplings from future experimental data. Results for h\ 2 couplings are similar to those 
obtained for the same values of h^ 4 . The only difference comes from the fact that the 
CP-violating couplings do not interfere with the SM. For a configuration like (b) with 
hj = hf = 0.01, hi = h% = 10~ 4 , the cross section is given by the average of both 
curves in the inset plot of Figure [5]. 

The situation is more complicated for production since the neutrino is not ob- 
served. Nevertheless, by identifying the transverse momentum of the neutrino with the 
missing transverse momentum, and assuming the W boson to be on shell, it is possi- 
ble to reconstruct the neutrino kinematics (particularly the longitudinal momentum) 
with a twofold ambiguity. In the case of the Tevatron, since it is a pp collider, it is 
possible to choose the 'correct' neutrino kinematics 73% of the times by selecting the 
maximum (minimum) of the two reconstructed values for the longitudinal momentum 
of the neutrino for W +r y {W~ r y) |26| . 

This is not true at the LHC, where due to the symmetry of the colliding beams 
both reconstructed kinematics have equal chances to be correct. Fortunately, in the 
case of anomalous couplings, we are interested in a efficient way to reconstruct the s 
rather than the full kinematics. Again there are two possible values of s. It turns out 
that there is a simple method to choose the 'correct' one 66% of the times at the LHC 
(73% of the times at Tevatron) by selecting the minimum s min of the two reconstructed 
values. This applies to both W +, j and W~^( production. Furthermore, we checked 
that the selected value s m i n differs in almost 90% of the events by less than 10% from 
the exact value s. This is likely to be precise enough, since the data will be collected 
in sizeable bins of s and the anomalous parameters are not expected to change very 
rapidly as a function of the energy. To quantify the advantage of the method, we 
show in Figures || and |7] the correlations of v^Wi with The cross section drops 
very rapidly for increasing y/§ — V^mm- This correlation also holds in the particularly 
interesting large y/§ region and also for the anomalous contribution. To investigate 
the latter point, we show in Figure [7] the correlation for (already experimentally ruled 
out) huge values of Ak = 0.8 and A = 0.2. For this figure we still use an ordinary form 
factor but in order to increase the anomalous contribution further we set A = 1 TeV. 

Finally, in Figure |8] we show the same correlation for production at the Tevatron 
(Run II), which corresponds to pp scattering at y/s = 2 TeV. We impose the same 
kinematical cuts as for the LHC, with the following exceptions: for the transverse 
momentum of the photon we require pj, > 10 GeV and the rapidities of the observed 
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Figure 6: The cross section for TF + 7 production (in pb/bin) as as function of 
V§ and Vsmin (in GeV) in order to illustrate the steep fall of a for increasing 
I — \/.§mm| • The left plot shows the full range of yi whereas the right plot 
concentrates on the particularly interesting large y/l region. 




Figure 7: The same as Figure || but with huge anomalous couplings Ak = 
0.8, A = 0.2 and A = 1 TeV. 
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lepton has to be in the range \rj\ < 1.5 (instead of 2.4). As for the LHC, there is a 
strong correlation between the 'true' center of mass energy \fi and the 'reconstructed' 
one v^mm- Again, this correlation also holds in the large region. In Figure [| no 
anomalous couplings are included, but we have checked that also for the Tevatron the 
inclusion of large anomalous couplings does not spoil the correlations. 





Figure 8: The same as Figure S but for Tevatron. 



As a result of these investigations we conclude that even in the case of produc- 
tion reliable bounds for anomalous couplings as a function of s can be obtained. To this 
end one merely has to do an analysis as in the Z7 case but with s m j n replacing s. This 
has the advantage that the anomalous effects can be quantized without introducing the 
ambiguity of form factors. Such a procedure would certainly facilitate a comparison 
of various bounds from different experiments. Finally, we note that in principle any 
quantity which has a very strong correlation with s can be used. However, we could 
not find any better candidate than s m i n . In particular, the correlations of s with the 



cluster mass and transverse mass respectively is not quite as strong [27 



5 Conclusions 



In this work, we have presented a general purpose Monte Carlo program for the calcu- 
lation of any infrared safe observable in Wj and Zj production at hadron colliders at 
next-to- leading order in a s . The leptonic decay of the W and Z-boson respectively has 
been included in the narrow-width approximation. We retained all spin information 
via decay angle correlations and thereby generalized previous calculations |9|, [10], pT| . 
We also included anomalous triple gauge boson couplings at NLO in a s and presented 
the analytical expressions for the corresponding amplitudes. 
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As an illustration of the usefulness of the program, we have studied several observ- 
ables for the LHC Generally we find that the NLO corrections are relevant for all of 



them, confirming results of refs. flq, || [LOL 11 



Moreover, we searched for the kinematical regions were the effect of anomalous cou- 
plings is amplified and proposed an alternative way to study its energy dependence. 
Using the strong correlations between the partonic center of mass energy and a mea- 
surable variable, s m j n , makes it possible to extract anomalous couplings from the data 
without need to introduce ad hoc form factors. Such an analysis is possible even for 
the Wj process and, in our view, should be undertaken at the Tevatron and LHC. 
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6 Appendix 

The helicity amplitudes that are needed for the calculation of and Z7 production 
at next-to- leading order can be found in ref. [0. In this appendix we list the additional 
amplitudes that are needed in the presence of anomalous couplings which lead to the 
non-standard vertices as given in eqs. (0) and (01). We use the notation and conventions 



of refs. [12 



We start with the W^f amplitudes. In order to maintain electromagnetic gauge 
invariance we set gj = and only allow anomalous couplings A/t 7 and A 7 . The only 
diagram that gets modified by these anomalous couplings is the diagram with a WW^ 
coupling. Fortunately, this diagram does not contribute to the rather complicated finite 
pieces F®,F^ of the amplitudes, given in eqs. (4.6) and (4.7) of ref. |i~2f . Therefore, 
only the tree-level amplitudes have to be computed with anomalous couplings. 

We do not have to list the explicit results for all possible helicities of the photon 
and the gluon. In order to reverse the helicities of all gauge bosons, i.e. the photon and 
the gluon (if the latter is present), we merely have to apply a 'flip' operation, defined 
as 

flip2 : 1 2; 3^4; (ab) «-> [ab] ; and A a <-> A b (for Wi) (6) 

The amplitudes A^' 3 ", A^ ee ' h given in eqs. (4.4) and (4.5) of ref. @ are modified 
as follows in the presence of anomalous couplings: 

a tree, a a tree, a , a5 

^ACri ~ ^7 ~+~ ^AC I ' I 

A tree, b /itree,b a 5 fo\ 

AA.C,7 — ^7 — ^AC \°) 
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where 

A\ c = ( ~ t[A5 \ n J (^ + A 7 ) (13) [25] [34] + A^(l|5|2) [45] ) (9) 
2s 34 (si2 - s 34 ) [34] v 



As usual, this will also lead to a modification of the one-loop amplitudes. The 
corresponding divergent pieces now read crVA^^ and crVA^Q^ respectively, where 



In the case of the bremsstrahlung amplitudes we have to consider two cases. The 
additional gluon can have positive or negative helicity. The corresponding Standard 
Model amplitudes are given in eqs. (4.9) to (4.12) of ref. [fTIJ. In the presence of 
anomalous couplings Ak and A they are modified as follows. 

Al^ 7 (l-, 2+ 3~ 4+ 5+ 6 + ) = A^(l-, 2+ 3", 4 + , 5+ 6 + ) + A? c (11) 

AK 7 (r , 2 + , 3~ 4+, 5-, 6 + ) = 4^(1- 2+ 3", 4+ 5", 6+) + Ajfe (12) 

At^(l-, 2+ 3", 4+ 5 + , 6 + ) = A^\r, 2+, 3", 4+ 5 + , 6+) - A 6 + c (13) 

AX5f 7 (l-, 2+, 3" 4+ 5-, 6 + ) = C b (l" 2+ , 3", 4+ 5" 6+) - Afc (14) 

where we defined 

A AC = » 9 2 /f ](l|2+ w!l/. 9 \ ( A7 ( 34 ) < 15 > ( 45 ) + ( A « 7 + ( 13 ) S 34) (15) 

2 S3W126 - (lo) (62) v > 

A e- s (35) (15) 

AC 2 S | 4 (t 126 - S3 4)(16)(62) 

(A 7 (35) (1|2 + 6|5> [34] - (A/^ + A^)(l|2 + 6|4)s 34 ) (16) 

We now turn to the amplitudes for Z7 production. In the Standard Model, there 
are no diagrams with a triple gauge boson coupling and the corresponding amplitudes 
^tree,s can s j m p}y Iqq obtained as the symmetric combination of the Wj amplitudes 

^tree,s ^tree,a _|_ ^jtree,b (17) 

(see eq. (4.15) of ref. fl^l). These combinations can be simplified somewhat but we 
refrain from listing the simplified versions. 

The amplitudes related to an anomalous ZZj or Z77 coupling will be denoted by 
^5/6 Ac - ^ n ^ ne f° rmer case, the intermediate vector boson is a Z-boson whereas in 
the latter case it is a 7. This results in different couplings of the intermediate vector 
boson to the initial state quarks. Apart from this difference, the amplitudes with an 
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intermediate Z and 7 are the same. The anomalous couplings always appear in the 
combination 

u z h u z h u z h u z h 

uZh - n i . h z h - ft 2 . iz/ 7 _ n 3 iz/ 7 _ n A 

hl ~ Mf ' h2 ~ Mf ' h * ~ Ml ' h4 ~ Mf ' [l *> 
We start with the tree-level amplitude for a positive helicity photon 

4% = 4^{( fh * h + ~ h " H ) 2 < 13 ) t 25 ] t 45 ] + (^? 7 + ^ /7 ) (12) (3|5|4) [25] 2 ) (19) 

In order to obtain the amplitudes for a negative helicity photon we have to apply 
the 'flip2' operation defined in eq. (|6|) in addition to 

~ h z h _~ h z h. ~ h z h ^ _~ h z h . (2Q) 

The anomalous bremsstrahlung amplitude with an additional positive helicity gluon 

reads 

= 4,34 (16) (26) X (21) 
({%hl h + h z z h )2 (13) [45] (1|2 + 6|5) + (fh z 2 h + /if /7 )(3|5|4)(l|2 + 6|5) 2 ) 
whereas for a negative helicity gluon 

A ^ = 4,34 [16] [26] X (22) 
({ihl h + h zh )2 [45] [25] (3|1 + 6|2) + {%h z 2 h + /if 7 )t 126 (3|5|4) [25] 2 ) 

Again, the operation 'flip2' reverses the helicities of the photon and the gluon. Finally 
we mention that in order to get the amplitudes with a positive helicity lepton, 3 + , we 
simply have to exchange 3 <-» 4 in the amplitudes presented above. Correspondingly, 
the amplitudes with opposite helicities for the partons are obtained by a simple 1^2 
crossing. 
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